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1 Introduction 



In this paper we continue our investigation of the thermodynamic Bethe ansatz (TBA) 
equations based on the quantum deformed AdSs x S 5 superstring S-matrix started in [1]. 
This story builds on the remarkably successful application of integrability to the AdS / CFT 
correspondence [2], relating the free AdSs x S 5 superstring to maximally supersymmetric 
Yang-Mills (SYM) theory in the planar limit [3, 4]. In this setting, the spectrum of scaling 
dimensions of gauge invariant operators in SYM is equal to the spectrum of the AdSs x S 5 
superstring - an integrable quantum field theory in finite volume. 

The finite size ground state energy of the AdSs x S 5 superstring can be determined 
from the thermodynamics of its mirror model [5, 6], resulting in a set of ground state 
TBA equations [7-9] based on the string hypothesis of [10]. 1 These equations imply the 
corresponding Y-system [13] combined with intricate analytic properties studied in [14, 15]. 
Excited states can be described through a type of analytic continuation of the ground state 
TBA equations [16-19], explicitly done for various states in [9, 20-24]. The testing ground 
for these ideas has been the Konishi state, where the TBA approach [25, 26] agrees with 
Liischer's perturbative treatment [27-30] and with the dual field theory up to five loops 
[31-33]. The TBA and Liischer based results were recently extended to six and seven loop 
order respectively [34, 35], awaiting gauge theory verification. 

The TBA approach has also been successfully applied to the cusp anomalous dimension 
in the context of Wilson loops [36, 37], and there have been interesting developments in 
describing the TBA through a finite set of non- linear integral equations [38, 39] 2 , see also 
[40]. In a specific limit, these last two topics merged nicely in [41]. More generic boundary 
conditions also appear to be within reach of the AdS 5 x S 5 TBA [42]. 

In our previous work, we derived ground state TBA equations based on the quantum 
deformed psu(2|2)® 2 invariant S-matrix [43, 44] when the deformation parameter q is a 
root of unity, more precisely q = e m l k with k an integer greater than two. This S-matrix 
is conjectured [44-46] to interpolate between the S-matrix of the light-cone gauge fixed 
AdSs x S 5 superstring at q = 1 and the S-matrix of the Pohlmeyer reduced AdSs x S 5 
superstring [47] 3 in a relativistic limit where the coupling g is taken to infinity. Because 
of the special nature of the deformation the TBA equations for this interpolating theory 
only contain a finite number of Y-functions with an interesting structure between them, as 
summarized schematically in figure 1. We will discuss the actual equations in some detail 
in later sections. 

We would like to extend these TBA equations to excited states in similar spirit to 
what we just described for the undeformed model. In the present paper we provide the 
necessary means to do so in the form of an asymptotic solution. This asymptotic solution 
is a solution of the excited state TBA equations in the limit where the momentum carrying 
Yq functions are small. Once found by for example the contour deformation trick [19], these 

1 For recent reviews in the present context see e.g. [11, 12]. 

2 In fact, the approach of [38] was used to compute the six-loop Konishi result of [34]. 
3 Recently the Pohlmeyer reduction procedure was also worked out for a string moving only in AdSs, 
which can be a more natural 'vacuum' depending on the physical context [48]. 
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Figure 1. The structure of the deformed TBA and Y-system. The green, teal (with ±) and blue 
dots indicate what we call Y M \ W , Y±, and Y M \ VW functions respectively; their coupling is nearest 
neighbour apart from the indicated coupling near the corner. The two wings of this diagram 
represent the two copies of the deformed Hubbard model that are coupled via the momentum 
carrying Yg-functions (purple) of the mirror theory. Note that Y\ couples to Y-, but no Yq couples 
to Y + in a local fashion. The lime-green bubbles signify the special state and twist dependent 



relation between Y u v \ w and Y t 



' k— 1| (v)w 



entering the Y-system. 



excited state TBA equations hold also outside this limit. Apart from the finite number of 
Y- functions and deformed integration kernels entering in the equations, this story mirrors 
the undeformed one perfectly up to this point. However, the fact that we have only finitely 
many Y-functions - a 'boundary' - introduces an important difference. 

In the undeformed case, we can derive a set of so-called Y-system relations between the 
Y-functions which appear to be of the same form for any excited state; they are universal 
[13]. Similarly, the undeformed Y-system is universal for all twisted boundary conditions 
in the string theory [49-53]. For our deformed model on the other hand, the Y-system is 
no longer universal in either sense. Near the boundaries of figure 1 the Y-system appears 
to pick up a mild dependence on (the excitation numbers of) the state under consideration 
as well as the twisted boundary conditions. For twisted boundary conditions this can 
already be observed for the XXZ spin chain [54], however we have not encountered this 
state dependence before. As this effect appears to depend crucially on the combination 
of the deformation with the fact that we have a nested system, perhaps this is not too 
surprising. 

As we will discuss in detail right after the introduction, this dependence comes about 
in two ways for either cause. Firstly the inverse relation between Y ^ w and Y k _i\^ w 
is not preserved by chemical potentials associated to twisted boundary conditions. This 
is because the 0\(v)w and k — l\(v)w strings have wildly different charges rather than 
the opposite ones which would be required to preserve the inverse relation. Secondly, 
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while twisted boundary conditions result in chemical potentials that are annihilated by the 
discrete Laplace operator used to derive the Y-system, near the boundary we necessarily 
apply a slightly different operator which leaves a chemical potential dependence. Next, 
the state dependence comes in firstly because the driving terms for Y u v \ w and Y k _iu v \ w 
naturally do not preserve their inverse relation; while resulting in kernels with opposite 
signs, their S-matrices are in fact anti-inverse rather than inverse. Furthermore, while the 
kernels are annihilated by the boundary operator we apply to obtain the Y-system, this 
only implies that the driving terms should be annihilated up to a constant term. This 
constant term turns out to be another minus sign in the game. These special relations 
between Yqu v \ w and Y k _i\^ w are emphasized by the lime-green bubbles in figure 1. This 
interesting dependence is of course reproduced by our asymptotic solution. 

The main results of this paper are the explicit derivation of the excited state Y-system 
and the independent derivation of its asymptotic solution, which showcases the highly 
involved structure of this interesting model while at the same time providing a nontrivial 
check of our string hypothesis. In addition to this we discuss the relativistic limit of 
our TBA equations, where they reduce to the TBA equations based on the conjectured 
S-matrix for the Pohlmeyer reduced theory. Of course in this limit the TBA equations 
simplify considerably. Most notably, the Y-functions become meromorphic functions on 
the u-plane and all nontrivial kernels including the dressing phase disappear completely 
from the simplified TBA equations. 4 

We have also carefully derived an explicit expression for the g-deformed mirror bound 
state dressing phase, and proved that it satisfies properties completely analogous to the 
undeformed case. Along the way we also show that the g-deformed transfer matrix used to 
construct the asymptotic solution can be obtained via fusion in the generating functional 
approach, and discuss the twisted ground state solution of the deformed model which takes 
a form very similar to the undeformed theory. 

This paper is organized as follows. In the next section we discuss the dependence 
of the Y-system on twists and excitation numbers, give a practical set of simplified TBA 
equations and the Y-system where these features have been explicitly taken into account, 
and discuss the mapping of the Y-system to the T-system in detail. Then in section 3 we 
derive the transfer matrix for the g-deformed theory with twists, which is used in section 4 
to construct the asymptotic solution of the TBA equations and Y-system. We then discuss 
the relativistic limit of our TBA equations in section 5, and finish up with a discussion of 
our results and outlook in section 6. Various background material and technical details 
are presented in appendices; in particular we give a complete derivation of the deformed 
mirror bound state dressing phase in appendix D. 

2 The Y-system for excited states and twists 

Typically the Y-system associated to a set of TBA equations is universal in the sense that 
it is the same for all excited states and independent of chemical potentials associated to 

4 The information contained in the kernels is of course still contained implicitly in the asymptotics of the 
Y-functions. 
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symmetries of the theory. This is not true for the present theory, a fact which is intimately 
related to having a finite number of TBA equations. 5 Let us begin by showing the effect of 
considering an excited state at the level of the canonical TBA equations, on the simplified 
TBA equations and associated Y-system. We will work this out explicitly for the simple 
case of w-strings. 

2.1 Excited states in the deformed mirror TBA 

The canonical ground state TBA equations for w strings are given by [1] 

log Y^ w = log ( 1 + i ) * K NM + log V * K M , (2-1) 



1 - 


1 
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1 - 


1 




'+ 



M\w & I 1 v (r) 

where we have not eliminated the Y- function corresponding to negative parity strings yet, 
and r = ± is a label referring to the two copies of su(2|2). Hence, it is important to note 
that the indices M and N are general and run over the length one to k — 1 positive parity 
strings, and the length one negative parity string. Denoting the negative parity length one 
string as type we note that the S-matrices the above kernels are derived from have the 
properties 6 

S k - 1 S° = -l, S Mk ~ 1 S M0 = l. (2.2) 

This immediately implies Y^ = for the ground state. For excited states 

however, we should be more careful. 

In general, the TBA equations for an excited state are the same as those for the 
ground state, up to a set of driving terms. These driving terms depend on the state under 
consideration and are different for each particle type's TBA equation. Let us for simplicity 
assume that we have a state with K w roots of 1 — Y_ that enter the TBA equations and 
no further roots. Considering the above form of the TBA equations, this would result in 
driving term contributions of the form 

5>gS w (rM-u-i), (2.3) 
i=i 

in the equation for \ogY^,{u). Because of the first relation in eq. (2.2), this results in 
a relative sign between the driving terms for logY^-ilw and loglo|uj> exactly in line with 
their Y-functions being inverse, as for the ground state. However, due to the minus sign in 
the first of eqs. (2.2) we also get an extra factor of K w m7 This means that while 

Y $ Y k-x\ w = 1 for even ^, (2.4) 



5 While atypical, chemical potentials (magnetic fields) already appear in the Y-system of the XXZ spin 
chain at roots of unity [54]. 

6 The explicit definition and properties of the kernels and S-matrices can be found in [1]. 

7 Of course the reader can argue that an S-matrix-like driving term log S(u,v) based on a kernel K is 
defined up functions of v only. This ambiguity should then be fixed by comparison to an asymptotic solution 
for the Y-functions. We do not have a direct asymptotic form for Yo\ w , making this impossible. However, let 
us firstly note that the chosen driving terms are the completely natural choice. Moreover, as we will discuss 
below the resulting equations are compatible with the available asymptotic form for the other Y-functions, 
while other potential forms of the driving terms are not. Hence there is no true ambiguity. 
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we have 

yPyP,, =-1 for odd K w . (2.5) 

\}\w k—l\w w <• ' 

Now let us come back to our original assumption. If we were to consider additional roots 
of Y± or 1 — Y+, clearly K w above should be replaced by the total number of all these 
roots, taking relative signs between driving terms into account. Furthermore, due to the 
plus sign in the second relation in eq.(2.2) additional roots of l + Y M \ w do not further affect 
this relation. 

Next, these driving terms also play a direct and important role in the simplified TBA 
equations and Y-system for Y k _x\ w . Conventionally, the TBA equations can be rewrit- 
ten in simplified form by noting identities between various kernels. As convolutions with 
convoluted kernels turn into convolutions with the standard kernel s, the driving terms 
in the simplified TBA equations turn into logS terms, where S is the standard S'-matrix 
associated to s 

Not surprisingly, when acting with s _1 to obtain the Y-system these driving terms are an- 
nihilated. 8 Up to potential subtleties with branch cuts, these identities between nontrivial 
S-matrices and the subsequent annihilation of the remainder by s _1 can be summarized 
by the statement that these S-matrices satisfy the discrete Laplace equation 

S M (u,v + i/g)S M (u,v-i/g) 

S M + 1 (u,v)S M - 1 (u,v) ' 1 J 

which holds for any of the relevant S-matrices in what follows. Since we have a boundary 
however, the identities we use there are equivalent to getting a final contribution of 

i ^^l^'-^ ^^M (2.8, 

in the simplified TBA equation for Y k _i\ w . Of course S k (u,v) is typically constant such 
that its associated kernel vanishes. However, with the auxiliary problem being 2k periodic, 
S k is not quite one, but rather minus one. This means that the simplified TBA equation 
for gets an extra factor of K w m jl. Repeating the above analysis for the canonical 

TBA equations for vw strings, 

1 _ _J_ 

log Y$ = log (l + * K NM + log -^*K M - log (1 + Y Q ) * Kg M , (2.9) 

V Y N\vw J 1 — 

it is easy to see that the relation between Y \ vw and Y k _i\ vw should be sensitive to both 
the roots associated to Y± and those to associated to Yq, as SxvSxv 1 = — 1 as well. Let 
us denote the number additional roots due to Yq by K vw . 

At this point we might wonder whether the number of these roots effectively follows 
a pattern. To understand and motivate this, we should realize that physically speaking 



3 s 1 is defined below in eqn. (2.14). 
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the difference between the driving terms for the 0|(v)u> and k — l\(v)w string comes about 
because of their different charges; they feel the presence of the excited state differently. 
This hints that the number of roots could be directly related to the charge of the state, in 
particular to their excitation numbers. 

To continue along this track, let us start by considering generic states in the exten- 
sively studied sl(2) sector of the undeformed theory. Based on various studies of the TBA 
equations for excited states in this theory [9, 20, 21, 55], it seems reasonable to assume 
that here the number of relevant roots is restricted to those of 1 + Yq corresponding to 
the fundamental excitations in the string theory. This would mean that K vw = K l and 
K w = in this sector. Next, if we consider states with nonzero with real momenta 
such as those studied in [22] the natural generalization of this statement is that K vw = K l 
and K w = -^jy)- 9 At the level of the number of roots this simple picture is not quite true 
anymore for for example an orbifolded magnon [50] or states with complex momenta in 
the su(2) sector [23] . However, we should note that the relation between the Y-functions 
only depends on whether K w and K w + K vw are even or odd, and that in this sense any 
currently investigated state matches the described pattern. Assuming this picture not to 
change qualitatively in our deformed model, we then have the natural conjecture that in 
general 



We might wonder whether or not the excitation number i£"f4 should make an appear- 
ance. To understand let us first look at the simple case of w strings again. In essence 
there the above story tells us that the 'length of the spin chain' associated to a particular 
particle type (nested level) determines the appearance of the minus sign. However, then 
we would immediately expect a minus sign for vw strings depending on whether K 1 is 
even or odd, apparently contradicting the above. The reason for the appearance of KV\ in 
this picture is that in deriving the 'canonical' TBA equations for vw strings the canonical 
TBA equations for w strings were used to exchange the presence of infinitely many Y M \ W 
functions for the presence of Y± [1, 10]. If this had not been done, the direct dependence 
would be on K l only, where the dependence on K F\ would come in when eliminating the 
Ym\ w functions from the simplified TBA or Y-system. Prom this point of view the presence 
of w roots (nonzero KV^) cannot result in any potential minus signs because there is no 
deeper nested level. 

While the precise form of the general relations (2.10) is a conjecture, it is completely 
natural to expect simple relations of this type by the physical picture that they are caused 
by the difference in charges of the respective string solutions. Of course, this conjecture 
is also supported by the asymptotic Y-system we will discuss in section 4. Finally, state 
dependence is clearly present by the above discussion, regardless of the precise form. 

9 Our excitation numbers refer to the st(2) grading. 



*0|w*fc-l|w 
0\vw^k— l\vw 
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(2.10) 
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2.2 Boundary effects for boundary conditions in the mirror trick 

As it turns out, the Y-system for our mirror model also depends explicitly on boundary 
conditions in the original model, which manifest themselves as chemical potentials in the 
TBA equations [5, 56]. Normally such chemical potentials drop out of the Y-system; here 
their effect does not completely cancel near the boundary of the Y-system. 

Incorporating nontrivial boundary conditions in the original model corresponds to 
insertion of a defect operator in the path integral for the mirror model, as illustrated in 
figure 2. This operator results in the appearance of chemical potentials to the right hand 




Figure 2. Quasi-periodic boundary conditions give a defect operator in the mirror theory. The 
quasi-periodic boundary condition denoted by the black discontinuity on the circle of the original 
theory (blue), turns into a discontinuity in the mirror time evolution (red), which is equivalent to 
insertion of a defect operator D in the partition function. 

side of the canonical TBA equations, worked out carefully in the present context in [52, 53], 
see also [51]. However, as the negative parity string has charge one in appropriate units 
compared to the charge k — 1 of the length k — 1 string, their chemical potentials are not 
quite opposite but rather sum to k units of chemical potential. This means the presence 
of a chemical potential results in a further modified relation between Y u v \ w and Y k _i^ v ^ w , 
namely 

Yo\ iv)w Y k . lKv)w = (-lfU +Kl )e- k » . (2.11) 

Moreover, when deriving the simplified TBA equation for Y k _i^ w we are left with a 
chemical potential contribution just as we were left with an S-matrix as in eq. (2.8). This 
results in a total factor of —}^{k[i + iirK l ^(+K 1 )) in the simplified TBA equation for 

Yk— l\(v)w 

Intuitively these boundary chemical potentials as well as the minus signs discussed 
above come about because while the 0\(v)w and k — l\(v)w strings scatter inversely in the 
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thermodynamic limit (as if they were each others anti-particles), they do not have opposite 
charges. 

2.3 The simplified TBA equations and Y-system 

The above discussion makes it clear that the simplified TBA equations and Y-system pick 
up an explicit dependence on the chemical potential and (the excitation numbers of) the 
state under consideration, in addition to conventional driving terms that generically still 
enter the simplified TBA equations for excited states. Part of this dependence comes 
about through the boundary associated with having finitely many TBA equations, while 
the other part comes in explicitly by eliminating Y u v ^ w from the TBA equations. We could 
choose not to eliminate Y w v \ w from the TBA equations, but since this would not remove 
the dependence on excitation numbers and chemical potential completely in any case, and 
we would have to keep extra careful track of the origin of driving terms, it seems best to 
eliminate Y ] u\ w . This way the TBA equations can be treated as in the undeformed case, 
with an asymptotic solution for the remaining Y-functions and the contour deformation 
trick applying directly and straightforwardly. The only modification is that we additionally 
have the manifest dependence on the excitation numbers and chemical potential. Let us 
note again that this dependence does derive properly from the contour deformation trick 
and the full set of canonical TBA equations as discussed above; it is simply easier to 
manifest it before continuing. 

With this in mind we can write down the simplified TBA equations for excited states 
with nonzero chemical potentials, up to standard log S terms following from the analytic 
properties of the independent Y-functions. 10 The full mirror model has two sets of auxiliary 
equations labeled by an index r = ±. Suppressing the index, for one copy these are given 
by 



logY M \ vw 


= log 


; (1 + Y M+l \ vw ){\ + lM-l|mo) * s - log (1 + Y M +l) 


* s 


+ S M ,i log 


is 

1-Y+ 


log 


= log 


; (1 + Y fc _3|,™,)(l + Y k _ Avw )(\ + e x ™Y k _ 1]vw ) * s - 


•l0| 


5(1 + Y fc _! 


)*8, 


loglfe_l| vlB 


= log 


I (1 + y k -2\vw) * S - log (1 + Y k ) * 8 - *p , 












l_y-i 






(2.12) 


logY M | w 


= log 


; (1 + Y M+1 \ W ){1 + Y M -i\ w ) *s + 5 M ,i log x 

1 Y + 


is , 




logYfc_ 2 |™ 


= log 


; (1 + Ifc-3| W )(1 + *fc-i|«,)(l + e Xw Yk-i\w) * s , 








logY fe _!| w 


= log 


;(i + n-2|«,)*a- *f, 









log Y ± = - » y - log (1 + Y Q )* Kf + log 1 * K M + log \ + Y ^ W * if^-i , 

1 + Y M]w 1 + Y k - 1{w 



For Y± we still write the canonical TBA equations, in principle they can be brought to so-called hybrid 
form [20] but we will not pursue this here. 
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where 



A.vw 



X 



(r) 



ku?) +iirK 



ii 

(r) 



The chemical potentials we consider derive from twisted boundary conditions in the original 
model. As in the undeformed model [53] they are naturally parametrized in terms of the 
twists a r and (3 r as summarized in table 1. 



A 


mV 


M\w 


2iMp r 


M\vw 


2iMa r 


y 


i(a r - P r ) 


Q 


iQ(a + + a-) 



Table 1. Chemical potentials corresponding to twists for the particle content of the deformed 
mirror TBA. 

The momentum carrying Yq functions couple two sets of the above equations through 
log Yi = log - * s + log Yl ^F)J -ks-Ms, 

log^Q = fog 



r=± 



Q — 1 1 vw 



(i + y _i)(i + y Q+ i 
iog^ = log r^g- ** + log II f 1 + ^) f 1 + 

x -» fc— I r=± ^ fc-i|ti«i/ \ k-i\vw / 



(2.13) 



Note that in the equation for Y& just above we might expect a chemical potential term 
because [Iq is not in the kernel of 5Qk — 2s5Qk-\. However, the resulting factor of 2m 
is exactly cancelled by the chemical potentials for a (+) and a (— ) l\vw string, cf. table 
1, whose canonical TBA equations are used to obtain the simplified equation for Y^ as 
discussed in appendix E of [1]. The Y-system can be immediately read off from these 
equations by applying s _1 to them 11 , where 



f os (u) = ]xm[f(u + i/g - ie) + f(u - i/g + it)] 



(2.14) 



Explicitly the Y-system is given by 



11 For Y± the story is more subtle. Y_ has a Y-system equation obtained via TBA identities, while Y+ 
has no local Y-system equation. 
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w strings 



^MT« = (i + ^| W )^ ir ^rJ (2-15) 

*M|w*Af|u; = + ^M-l|w)(l + *Af+l|u;) > (2.16) 

^-21^^-21^ = C 1 + n-3|«,)(l + + e Xu, n-i|») , (2-17) 

Cii«n-n« = e " x ™ (! + n-2i«) • (2.i8) 



vw strings 



1 + Y 2\vw ( 1 - Y- \ e,(uB_|u|) 
y i>«, y i|^ = J + y 2 \1-Y+ J ^ 2 ' 19 ^ 



J- T -12 \ ± — J + / 
_ + _ (1 + ^M-1|dw)(1 + ^M+l| 

fc— 3|i)to) (1 

1 + *fc-2|mu 



Af|trcu M|™ j + Y M+1 

Y + Y~ = - " ' 

k—2\vw k— 2\vw ^ _|_ ^ 



1 + Ym+i 

(1 + Y k _ 3{vw )(l + r fe -i|^)(l + e^-n 

Y + y- - r -xvw 1 + Y k~ 2 \™ ( 2 no) 

1 k-l\vw x k-l\vw ~ e 1+lfe ' ^.^J 



y strings 
Q particles 



- - " i + ni„ i + i'i ■ ( ' 



n fi+^ 



Y+Y- Y {r) 
*Q*Q _ r =±\ / Q-l|t>«>- 



(2.25) 



Y Q+1 Y Q ^ (1 + Y Q - 1 )(1 + Y Q+1 )' 

(2.26) 



Y+Y- n ± ( 1 + ^)( 1 + ^ 

*k I k r— ± \ k—l\vw' \ k— l\vw ■* 



Y li 1 + Y k-i 

The above equations have an explicit dependence on \ that cannot be removed by 
simply redefining the Y-functions. However, let us note that by the redefinition 

X(v)w 

Y k-i\(v)w ~> e ~ 2 Y k-i\(v)w (2.27) 

the above equations do become manifestly symmetric in %, and in particular in the twist. 
Since we would expect the spectrum to be symmetric in the twist, this is only natural. 
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Figure 3. The T-hook lattice for T-functions satisfying the Hirota equations. There is a T- 
function T a ^ s indicated by a gray dot at any lattice site (s,a). The bigger coloured dots cover 
underlying gray dots and signify the Y-functions constructed out of nearest neighbour T-functions; 
the stars denote the special Y-functions which have no direct construction in terms of standard 
T-functions. By eliminating the special Y-functions, the consistency of the Y-system near the 
boundary is equivalent to a set of non-standard additional relations between the T-functions in 
each respective yellow bubble. These relations allow us to consider a truncation of the in principle 
infinite lattice of T-functions satisfying the Hirota equations. 



2.4 T-system 

It is useful to express the Y-functions in terms of so-called T-functions. In the undeformed 
case the Y-system equations then become equivalent to a set of Hirota equations for these 
T-functions. Here we will study what happens in the quantum deformed model. 

As in the undeformed case we relate our Y-system to the T-system by making the 
following identifications [13, 20] 

v (±) _ t i,±(m+2)Ti,±m (±) _ T m +2,±iTm,±i v _ T M ,iT M -i , v 

M\W T T 1 Y M\VW rji rp ! 1 M — rp rp 

y(±) _ T2 }±1 T2,±3 ^ y (±) _ T 0>±1 T 2)±1 ^ ^ 2g ^ 

+ 7l,±27l,±3 ' ^1 j oTi ) ±2 

In the undeformed limit, the Y-system equations are equivalent to the Hirota equations 

Ta, s Ta )S = T a +l,sT a -l,s + T a .s+lT a ,s-l- (2.30) 

plus the boundary conditions that T as vanishes if both a and \s\ are greater than two. These 
boundary conditions mean that the T-functions lie on a so-called fat hook, as illustrated in 
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figure 3. Away from M = k — l,k — 2 this identification and equivalence also immediately 
work for our g-deformed model. However, at the boundary where a,s = k, k+1 we obviously 
find special relations. In principle we can eliminate the T-functions that do not satisfy the 
standard Hirota equations to obtain a set of standard T-functions that solve the Hirota 
equations 12 , plus a set of complicated additional relations signifying the closure of Y-system 
at M = k — l,k — 2. This point of view is illustrated in figure 3, where these extra relations 
are schematically indicated by the yellow shading, simultaneously marking the boundary 
of our T-system. Note that these relations are not between nearest neighbours only, as the 
Hirota equations are. 

Before eliminating the 'non-Hirota' T-functions, the Y-system is equivalent to the bulk 
Hirota equations plus additional equations for Ti^r^u, T&+i,±i, and T^+ifl- Firstly, it is 
easy to see that the Y-system equations for w strings for M = k — 2 and M = k — 1 give 



^lfc^l k ~ (^o,fe^2,fc + 7l,fc-l2i,fc+l) 



1 + e 



+i 



rp+ rp- 

-nJM-l l,fc+l 



.,(+) 



e~ x - T 0Ml T 2 , k +i 



(2.31) 
(2.32) 



Note the dependence on the chemical potential, and that the first factor on the left hand 
side of the first equation is of the usual Hirota form. Similarly, the last two equations for 
■uu>-strings are 



,(+). 



T, 



fc+i,i fc+i,i 



Xvw 



-./+) Tu_ i iTi,t | i 

1 _j_ gXvw 1 ' 1 I 



Tk+i,oTk+i,2 ■ 



(2.33) 
(2.34) 



Of course there are similar equations for the (— ) wing. In addition, the equation for 
gives 



rp+ rp- 



T^2 rp2 

± k+l,l ± k+l,- 



2,0 



1,0 



k.2 



'-k,-2 



[ k,0 



(2.35) 



These equations together with the bulk Hirota equations define our T-functions on a finite 
- and again truly fat - hook, as illustrated in figure 4. Restricted to the auxiliary problem 
the functional relations satisfied by the T-functions only involve the conventional nearest 
neighbours. The equation for Tfc+i o necessarily involves non nearest neighbour terms 
however, and as T^+ifl already enters in a regular Hirota equation it does not appear to 
make much sense to attempt to change this by introducing further auxiliary T-functions. 
The parametrization of Y-functions in terms of T-functions is left invariant under 



T, 



a+s] [a—s] [— a+s] 



a, s -> g[ g 2 % 



9a 



-L a 



(2.36) 



where the gi are arbitrary functions, presenting a gauge transformation on the T-functions. 
Because the Y-system is invariant, the T-system is invariant under these transformations as 



12 Since we have only a finite number of T-functions by construction, here we mean that they solve the 
Hirota equations in the 'forward' sense, going outward from the origin of the lattice. 
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(0,0) 



Figure 4. The finite T-hook lattice for T-functions. Also here there is a T-function T QiS indicated 
by a gray dot at any lattice site (s, a). The gray triangles indicate the 'non-Hirota type' T-functions. 
By keeping these T-functions the nontrivial boundary consistency conditions are equivalent to the 
set of relatively simple nearest neighbour equations (2.31)-(2.34). These relations are indicated 
by the wavy lines, while straight lines indicate the standard Hirota relations. Equation (2.35) for 
Tfe+i.o remains non nearest neighbour, involving the T-functions in the yellow bubbles. 

well. This in particular includes the special equations (2.31)-(2.35) as can also be verified 
by direct substitution. 

Asymptotically the 'Hirota type' T-functions coincide with the transfer matrix based 
on the psu g (2|2)® 2 invariant S-matrix. It should be clear however that Ti^^^i^fc+i^l 
require some care; we will see this explicitly below. Let us now derive the transfer matrix 
and then discuss how it provides us with the asymptotic solution to the T-system. 

3 Transfer matrix and fusion 

In this section we derive the (twisted) bound state transfer matrix from first principles. 
We will also show that it can be obtained by a fusion procedure. We will work on the 
u-plane (additive parameter). In this presentation the parameters describing a M- 
particle bound state can be conveniently defined via x ± = x(u± iJ f). More generally, let 
us introduce the notation x^ = x(u + y), then these parameters satisfy 

^ w+ J0-4 M+ ^)4°-?)H)' <"> 
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where 

(3.2) 



5 sin' k 
1 + g 2 sin 2 £ ' 



The derivation presented here will be very similar to the ones presented in [53, 57] and we 
will mainly restrict ourselves to highlighting the differences. 

3.1 Algebraic Bethe Ansatz 

General procedure. Let us first briefly outline the general procedure. We consider K l 
fundamental particles with rapidities ui, . . . , u K i. To these particles we add an auxiliary 
a-particle bound state with rapidity u. This system is described by the following tensor 
product space 

V := V a (u) ® Vi(ui) ® . . . (g> V K i(u K i), (3.3) 

V v ' 

Vp 

where Vi is the fundamental representation with rapidity U{. We divide the states into 
an auxiliary piece V a and a physical piece V p . The monodromy matrix is then defined as 
follows 

K 1 

T a (u, {m}) := Y[&ai(u,Ui). (3.4) 

1=1 

This operator manifestly depends on the representation of the auxiliary particle. 

The monodromy matrix is a matrix in the auxiliary space V a {u) with matrix entries 
being operators acting on Vp. Taking a basis |e/) (see (A. 15)) for V a (u) and a basis 1/^) 
for Vp, the action of the monodromy matrix T = Ta(u, {u^) on the total space V can be 
written as 

T(\ej) |/a)) = E T /a M ® \f B ). (3.5) 

J,B 

The matrix entries of the monodromy matrix can then accordingly be denoted as 

r\ ei ) = J2T/\ej) , (3.6) 
J 

while the action of the matrix elements Tj as operators on Vp can easily be read off to be 

r/|/ A ) = ^r//|/ B ). (3.7) 

B 

The operators Tf satisfy non-trivial commutation relations among themselves. Consider 
two different auxiliary spaces V a (u a ),Vb(ub). The Yang-Baxter equation for the S-matrix 
implies that 

$ab(u a , u b ) T a (u a , {Ui})Tb(u b , {ui}) = T b (u b , {ui})T a (u a , {ui}) §> a b(u a , u b ) , (3.8) 
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where §> ab (u a ,u b ) is the S-matrix describing the scattering between the two auxiliary par- 
ticles. By explicitly working out these relations, we can find the commutation relations 
between the different matrix elements of the monodromy matrix. These fundamental com- 
mutation relations (3.8) constitute the cornerstone of the Algebraic Bethe Ansatz [58]. 
The transfer matrix is subsequently defined as 

t a (u, {u t }) := sti a T a (u, { Ui }) = (-1) W T/, (3.9) 

and it can be viewed as an operator acting on the physical space Vp. By virtue of the 
Yang-Baxter equation it defines an infinite set of commuting quantities 

[ta(Ua),t b (u b )}=0, (3.10) 

making the integrability of the model explicit. Our goal is to find the eigenvalues of the 
transfer matrix by means of the algebraic Bethe Ansatz. 

Twists. Since the coproduct of the Cartan generators Hi is unaffected by quantum trans- 
formations, the group-like element G = exp^iaHi+ifSH^) gives a symmetry of the S-matrix 
in the sense that 

[G a ®G b ,S ab ] = 0. (3.11) 

We can then introduce a family of twisted transfer matrices 

t°(u, { Ul }) := str a G a T a {u, {m}). (3.12) 

As a consequence of (3.11), the twisted transfer matrix preserves the integrability property 
[ta(u a ),tf(u b )] = 0. Such a twisted transfer matrix describes an integrable model with 
quasi-periodic boundary conditions [59]. 

3.2 Transfer Matrix on the Vacuum 

The algebraic Bethe ansatz basically consists of two ingredients. One defines a vacuum state 
and a set of creation operators that are elements of the monodromy matrix. Eigenstates of 
the transfer matrix are then build by acting with the creation operators on the vacuum. Like 
always, the transfer matrix respects the two su q (2) symmetries of the symmetry algebra 
of the S-matrix. In particular, the transfer matrix commutes with the Cartan generators 
iii,3 and consequently, the eigenspaces of the transfer matrix are labelled by the quantum 
labels K , K 111 corresponding to the eigenvalues of H\ 3. Explicitly, for a state in V p of the 
form 13 

K l 

§Z)\mi,ni,ki,k), (3.13) 



we have 



K l K 1 

K u = ^2mi + ni + 2ki, K m = ^ m + h. (3.14) 



8=1 1=1 



3 See appendix A.l for details and notation regarding the fundamental representation. 
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These numbers once again count the total number of fermions (K 11 ) and the total number 
of fermions of a specific species (K 111 ). 

In our case, we define the vacuum to be the highest weight state corresponding to the 
quantum numbers K 11 = K 111 = 

K 1 

|0) =®|0, 0,0,1). (3.15) 

i=i 

From the explicit form of the monodromy matrix (3.4) it is then straightforwardly seen 
that this is an eigenstate of the transfer matrix, with corresponding eigenvalue written in 
terms of S-matrix entries 

K 1 K 1 4 a-1 ia(a _2m-l) & 

A =e™ J] ^\u, Ul ) + J] (u, u t ) - £ £ ]J u,)+ 

i=l i=l 7=3 m=0 i=l 

a-1 

+ E E e " (a ~ 2m) ^^«,«l^^ (3-16) 

ai=l,3 m=l 

The last term can be mapped to a 2 x 2 matrix eigenvalue problem just like in the unde- 
formed case [57] involving the eigenvalues 

A _ x+ - xj x- - (x^)^ sinh %{u -u l + i^Dl) 
± i+-i:i+-( i M'»])±i sinh M( u _ Ui+ 

From the explicit S-matrix entries listed in Appendix B.l it is then straightforward to 
arrive at the following the eigenvalue of the transfer matrix on the vacuum 

K l _ _ j + a-1 / K l K l \ 

= e ^ a + e -^ a n g a X + 1 x - t _ x + + e e ia(a ~ 2m) n A ++n A ( 3 - 18 ) 

i=l X X i x+ X i m =l \i=l i=l / 

a— 1 if 1 2m-i + + «mh ""ff f,/ „.. I (a-2m-l)i x 

_ ( e tf + e -tf) ^ e «Mm-l) Jj g 2 x+-x^mh^(n-^+ 
m=o i=i J*[ x+ - X i smh^§(u-u; + 

In the g — ?• 1 limit it immediately reduces to the result from [53]. The untwisted transfer 
matrix simply follows from setting a,/3 = 0. 

3.3 Creation operators 

General states are constructed with the help of creation operators. Consider a monodromy 
matrix where the auxiliary space is in the fundamental representation. Let us denote the 
basis vectors as \ei), with i = 1,2,3,4. Then we define the following creation operators; 
two fermionic T^,T^ and one bosonic T\ operator. They have a clear interpretation. By 
acting with on the vacuum state you create an excitation of type a. More precisely, the 
creation operators increase the quantum numbers of the state as listed in Table 3.3. 

Notice that any state with quantum numbers K 11 , K 111 can be constructed by apply- 
ing these generators. Moreover, we see that both T3T4 and T\ have the same quantum 
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AK U 


AK m 


rpl 


1 


h- 1 


rpl 


1 





rpl 


2 


1 



Table 2. The quantum numbers of the creation operators. 



numbers, which indicates that they will mix. This has its origin in the presence of the 
scattering process that mixes bosons and fermions. In other words, a generic state will 
be a linear combination of all the operators that create the same quantum numbers. Of 
course, this can only be an eigenstate for specific values of the relative coefficients. 

The last crucial ingredient that is needed are the commutation relations. By computing 
the commutation relation between two creation operators, a nested structure is revealed. 
The fermionic creation operators can be commuted via an su(2) invariant S-matrix, r, in 
the following way 

Tl{ Ul )Tl(u 2 ) = - X \~ X { U T ^l rl^u, - n 2 )Ti(n 2 )T 7 1 (n 1 )+ (3.19) 

+ % {xt " 1 " r)( + 4 _ " x ' 2) U ^{n(u 2 )Tl( Ul ) - r.Vi^M), 

q e a /3 l _ x Jx 2 7172 

where a, p = 3, 4 and we take 634 = 1. The non-zero components of r(u) are rg| = r|| = 1 
and 



34 _ 43 _ sinh 2F „34 _ sin f e-^sin 



43 -' 34 -sinhS(.-f)' ' 34 -sinhS(.-f)' ' 43 -sinhf(u-f; 



(3.20) 



and it is simply the S-matrix of an inhomogeneous XXZ spin chain. As such it is unitary 
and satisfies the Yang-Baxter equation. 

The relevant parts of the other commutation relations that deal with the diagonal 
components of the monodromy matrix, and hence are important for the transfer matrix, 
are listed in Appendix B.2. 

3.4 Full Transfer Matrix 

In [57] the eigenvalues of the transfer matrix were computed quite generally in terms of 
elements of the S-matrix in (3.42). Indeed, it applies here as well, but the auxiliary problem 
is now described by (3.20) rather than the XXX spin chain. Let us introduce a short-hand 
notation for the different types of functions that appear 



r 9 ' n R ± _ TT 1 + 7 ] 

i£i / EF ' A J_ J. *£k r= 

n=l y 4 V« n =l Q 4 V x n 



OI.V- -l ' — /1 '-A 11 A£l ' ~A 11 A^l /—= > 

— " H ft v \/i?x . n 1 \ I Xn 

(3.21) 
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where u- 11 ^ = Vi and uf 1 ^ = Wi. We then arrive at 



^ _ jgg Va TT | V% ~ x / % I \ ~* Jl n -2m J a-2in , \ ^ rt g-2m- P q -2m 



p+ R~ i i 7.. V rj.— I / 021001 Z—s e 1iam /Oil 

-"•a -Da i=1 y« x v X ^ m= o e V a -2m m=1 e V a -2m 

0-1 p- p- 

»„-2,;, i -'.,-2,„ 



Va-2m ^a-2m-2 



... =0 V a -2m-l 

while the auxiliary parameters y m and w; m satisfy the auxiliary Bethe equations 



(3.22) 



1=e ^n^/|n :;:s;ii_:; + ij ; ^ 

f=i sinh 1 (™« - «i - \ ) ^ sinh f§ (w n -Wj + f)' 

Vj = ~ "I 1 (3.25) 

According to formalism of the analytic Bethe Ansatz these follow from the pole structure 
of the transfer matrix. The main Bethe equations pick up a a twist factor e ia as in 



where 



3.5 Generating Functional 

There is an alternative way to derive the transfer matrix in higher representations from 
the fundamental one. This method goes under the name of generating functional [60-62]. 

Up to an overall normalization, we can write the fundamental transfer matrix as a sum 
of four terms 

h l = e ™ + e - ia R -i B ~iQi _ etf^TQg _ e -V R l Q ^ Ql} l = t x + t 2 - t 3 - U- (3-26) 

Consider the shift operator D which is defined via Df(u) = f(u — ~)D, then we can 
introduce the generating functional W 

H ' = T3lkl5 [1 - 044D1[1 - 0t3D1 ra^ (3 ' 27) 

The above defined functional W then generates the different transfer matrices by expanding 
in the shift operator D according to 

W = Y,D a t a ^D\ W- 1 =J2 DS h,sD s . (3.28) 

a s 

The fundamental transfer matrix agrees by construction. It is then straightforward, but 
tedious, to match the transfer matrix term by term. 

Let us also spell out the transfer matrix t\ tS , which can be computed via the inverse 
generating functional according to (3.28). This transfer matrix will prove to be useful in 
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defining some Y- functions later. The transfer matrix is generically denned up to normal- 
ization, so let us adapt a convenient one for this purpose by setting 



*3 



(3.29) 



Then it is straightforward to work out W 1 explicitly 



o m 



-iQ 



C>+ R- s-1 p 2i/3mfrll frill frill 
XT ' ■•' — 2m J. — tf > — I 



III 

s-2m+l 



o+ * j/3(2m-lWI -QUI -QUI 



m=l 



Q 



-. -2m-lQ l s-2m+lQs+ ! 



(3.30) 



- C+ 



in 



III /-) ln 

s-2m+lV s +l 



m=0 



frill 



s-2m+l 



where we used that Rt s B^ s = Rj Bj . A different, related procedure is that of dualizing 
the transfer matrix. Details are provided in Appendix C where we compare ti iS to the 
dualized version of t a i . 



4 Asymptotic solution 

Now that we have derived the twisted transfer matrix, we can construct a solution for 
the Y-system when Yq —> Yq ~ 0. The Y-system in this limit is called the asymptotic 
Y-system, as it arises for example when J is large and the asymptotic Bethe Ansatz should 
describe the problem exactly. In this section all Y-functions will be understood in the 
asymptotic limit and to avoid cluttered notation, we will simply identify all Y-functions 
with their asymptotic limit and drop the superscript o. 

4.1 Asymptotic Y-system 

Let us first summarize the asymptotic Y-system equations for the different types of strings. 
As before, we only label the left and right copies of su g (2|2) explicitly in the equations for 
Q-particles. 

w strings. The equations for w strings do not contain the Yq functions so their form is 
simply unaffected 

(\ _ Y -1\ 9( - UB -\ U \) 

Y it Y i\ w = (! + y 2|J [tZ^) (41) 

Ym\ w y m\v, = ( 1 + y m+i\w)0- + Y m . 1]v> ) , M = 2,...,k-3, (42) 

Y k + -2\ w Y k Z 2lw = (1 + n-3| J(l + Y k _ 1]w )(l + e^Y k _ llw ) , (4.3) 

Y k t llw Y k -_ llw = e-^(l + Y k ^ w ). (44) 
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vw strings. By dropping the Yq terms we get a similar set of equations for uw-strings 

Y iU Y i\™ = a + Yq™) (t^y;) (4 - 5) 

Y M\vw Y M\vw = + Y M+X\vw)Q + ^M-l^w) > ( 4 -6) 



M\vw M\vw 

K-2\vw Y k-2\vw = ( 1 + Y k-3\vw)(l + Y k-l\uw)(l + e X ™ Y k-l\vw) , (4-7) 

Y kU™ Y k-n™ = e ~ Xvw{ ± + Yk -^ ■ ^ 

Apart from the first equation they are exactly of the same form as those for w strings. 
y strings. For y particles there is only one local equation. Asymptotically it reduces to 

Y + Y- = 1 + Yl \™ . (4.9) 

Q particles. The Y-system for Q-particles involves ratios of Yq functions. These survive 
in the asymptotic limit, resulting in the following set of equations. For Q = 1 we have 



!! IF 1 

r=± 

Then for Q = 2, . . . , k — 1 we find 



i?-=n i 1 -^)) ' f° r \ u \ < u b ■ (4.io) 



y+y- 
Q Q 



V V --II 1 + -(r)— ' (4U) 



-l\vw 



while at the boundary Q = k 



Y+Yr T-r / 1 \ / e-*«* 



n i+^r- ■ <« 2 > 



A; A; 

"yj" "11 F (r) IF v (r) 

fc-l r =± y z fc— l|t«a / \ 1 k-l\vw / 

4.2 Solution via transfer matrices 

In the asymptotic regime we can find a solution of the Y-system in terms of transfer 
matrices. The T-system indicates how this is done. We started out with a model with 
two copies of su g (2|2) symmetry. These copies form the left and right wing of the fat-hook 
and this structure should be manifest asymptotically. Hence, we are able to express the 
T-functions in the left and right hook in terms of the su g (2|2) transfer matrices that were 
derived in the previous section. Of course, the left and right copies can in general have 
different twists a±, (3±. 

The T-functions are defined up to gauge transformations (2.36), which we can use to 
set Tq^ = 1. Then, for a < k we make the following identification 

T a ,±l -> Ta,,xta,l\a=a±,P=P± ' T afi ~^ ( 4 - 13 ) 
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where T a ,i is simply a normalization of the transfer matrix such that 



8=1 

where the itj are the asymptotic rapidities of the excitations under consideration, 

Sj(2) = (S gi r 1 (E gU )~ a , ( 4 -!5) 

and S^ 1 * is the improved mirror bound state dressing phase with its second leg continued 
to the 'string' line; see appendix D for an explicit expression for and a discussion of the 
properties of E. The reason for this factor is that the transfer matrix (3.22) is constructed 
from the canonically normalized S-matrix; the S S U2) part of Tq simply reintroduces the 
correct dressing for the S-matrix. Secondly we would like to point out that e~ J£ -Q satisfies 
the discrete Laplace equation, so it is basically a gauge transformation. We see that this 
factor implies 

Y Q ~e- M Q, (4.16) 

which makes explicit that Yq is small for large J. Moreover, we now find that the asymp- 
totic Bethe equations are equivalent to 

Y u ( Ui ) = -l, (4.17) 

as in the undeformed model. 

Bazhanov-Reshetikhin formula. In section 3 we have explicitly constructed the trans- 
fer matrix t a ^i for short representations corresponding to a Young diagram with a rows and 
one column. However, (2.28) also contains T-functions with different indices. There is a 
well-known determinant formula that provides a formal solution of the Hirota equation in 
terms of t a \ as 

t a s (u)= det t a+m - n l (u + i(s + 1 - m - n)/g). (4.18) 

Km,n<s 



Equivalently, we could also use t\ s to generate all required transfer matrices 



14 



Asymptotic solution. It is now straightforward to derive the solution of the asymptotic 
Y-system in terms of transfer matrices. We get the following identification for M < k — 1 
and Q = 1 , . . . , k 

Y M\w — — -. > Y M\vw — — > Y Q — ±QtQ,ltQ-i 

J2.M+1 tM+1,2, 

Y_ = ^i, Y + = - t ^. (4.19) 

h,2 *3,2*1,2 



14 We remind the reader that ti >s in eqn. (3.30) is derived with a differently normalization from t a ,i in 
eqn. (3.22). The Bazhanov-Reshetikhin formula allows us to express ti. s in terms of t a ,i and doing so we 
find agreement with (3.30) up to this normalization factor. 



- 22 - 



As always, the auxiliary space uses the mirror parameterization, while the physical particles 
use the 'string' variables. 

However, this standard construction clearly cannot hold at M = k — 1 since we know 
that Ti^ + i and Tfc+i 1 satisfy a special set of equations. To correctly take the contributions 
of these functions into account, let us set 15 



Yi 



k— l\w 



Y, 



k— l\vw 



tk-1,1 



where 

tl,s 



Qi 



in 



■ R + 



8+1 

2 E 



and 

t a ,i = ( 



r+b: 



ia(3a+2) _Qa_ 
Ra B a 



8+1 

TP E 

na m=s+l 



e i/?(2m-i)gn_ 2m Qin_ a gni i 



s-2m-l^s-2m+ 



Ql-sQs-1 



iQ 



in 



+ 



(4.20) 



(4.21) 



s— 2m- 



in n 111 

s-2m+lV s +l 



a m= 



s— 2m- 



-iQ 



in 

s-2m+l 



X 11 
i=l 



m=a+l 



E 



R-a-2mB a -2m 



4 e ia(2m+l) Qll 



ill 



-2m-l 




R-a-2mB a -2m 



2m 



E 



^a-2m^a-2m 



2m+l , -i/3^a-2m-3 
I 6 



in 



V a -2m 

(4.22) 



Qa-2m Qa-2m-2 

Both terms contain only four terms and as such appear to be some special case of the 
fundamental transfer matrix. They are obtained by extending the sums that appear in the 
general expressions (3.22) and (3.30). 

From the explicit expressions for the transfer matrix, it is straightforward to check that 
(4.19) and (4.20) indeed provide a solution to the asymptotic Y-system. Let us highlight 
two reasons why this is quite non-trivial. 

First, the function Y^_iu v \ w a priori satisfies two different equations which relate t to 
t, and it is not obvious that both equations are compatible. Their compatibility crucially 
depends on the fact that q k = —1. 

The second non-trivial consistency check lies in the equation for From the asymp- 
totic equation for Yjj_i we can immediately deduce that Y^ is in fact given by the bulk 
expression at Q = k, while it satisfies an equation of a different type itself. The fact that 
equation (4.12) is satisfied nonetheless relies crucially on the factor Tq. Indeed, since Tq 
satsifies the discrete Laplace equation 



Q Q 



0+1 l£ 



(4.23) 



we see that in the RHS of (4.12) the dependence of Tq does not drop out. Following the 
arguments from Appendix E from [1] we can derive the following identity 



n-i 



(a 



10\2 
xyl 



K 1 

IL 

1=1 



x [k+l] 
x \k+l] 



1 



x [k 1] 3 
1 

x \k— 1] 3 



(4.24) 



This non-trivial factor makes (4.12) hold. 



5 In (4.20) ti, s is given by (3.30). 
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4.3 The twisted ground state solution 

As shown in [1], the ground state TBA equations of the deformed model have almost the 
same solution as the undeformed model [63]; only the boundary Y- functions are different. 
Interestingly enough it turns out that this feature persists for the twisted ground state. 
In the asymptotic limit, the constant bulk Y-system for the auxiliary problem is given 

by 

Yl = Y 2 = , (4.25) 

y m\ w = (1 + Y M -i\ w )(l + Ym+i\w) , (4-26) 
Ym\vw = (1 + Ym-i\vv>)(1 + Ym+i\ vw ) > (4.27) 

where potential M = terms on the right hand side are zero. The general solution to these 
recurrence relations for the bulk Y M \ W and Y M \ VW is given by 

Y M[w = [M] elb [M + 2U, (4.28) 
Y M \ VW = [M] eia [M + 2} e , a , (4.29) 

where a and b are undetermined complex numbers for the moment. Of course, we have yet 
to take into account the nontrivial boundary of the Y-system. These boundary equations 
read 

([k - 2} elb [k] etb ) 2 = (1 + [k - 3] eib [k - l] eib )(l + Y k _ 1]w )(l + e k ^Y k _ 1]w ) , (4.30) 
Yk-l\ w = e~ k ^(l + [k- 2] eib [k] eii ) , (4.31) 

with a similar equation for Y k _i\ VV} . These equations do not have a unique solution for 
arbitrary H\u v \ w and b(a). However, insisting that these equations satisfy the canonical 
TBA equations we immediately find that a and b are precisely a and (3 respectively as in 
table 1, and we get 

Y k _ 1]w = e- ikfS [k-lU, (4.32) 



( k-l\w 

Y k ^ llvw = e- ika [k-l} eia . (4.33) 



With this solution we know Y± by the Y-system equation above and can then find the 
Yq functions by their canonical TBA equations. The resulting Yq functions are formally 
identical to the Yq functions of the undeformed twisted ground state and are given by 

Yq = ([2] e , a+ - [2] e «, + )([2] eto _ - [2] e ^)[Q]^ + lQb«-e- JiQ{p) , (4.34) 

As the entire solution except the boundary w and vw Y-functions formally coincides with 
the undeformed twisted ground state solution [52, 53], this solution of the deformed TBA 
equations smoothly turns into the twisted ground state solution of the undeformed model 
in the simplest way possible. Computing the ground state energy for these solutions at 
leading order in g gives complicated expressions; instead of giving them we have plotted 
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E 




Figure 5. The ground state energy of the twisted deformed theory. The plot shows the energy 
as a function of k, for J = 2,3,4 with arbitrarily chosen twists a+ = ^j, a_ = /3+ = jjj, and 
/3_ = — 4j- The energies at fc = 1 and k — 2 are the continuation of our results valid for k > 2. 



the energies for some generic parameter values in figure 5. The energy converges to the 
undeformed result as p at large k. Explicitly 

E(k, J) = E{J) - ^ (E{J) - ^j^fiJ ~ 1)) + O (p) , (4.35) 
where E(J) is the ground state energy of the undeformed twisted model [53]. 



5 The relativistic limit 

Keeping k fixed and taking g — > oo in an appropriate fashion, we obtain a set of TBA 
equations based on the conjectured S-matrix of the Pohlmeyer reduced superstring, a semi- 
symmetric space sine-Gordon (ssssG) theory. The appropriate limit to take is to rescale 
our rapidities u — )■ | and take g — > oo keeping u fixed. In what follows by conventional 
abuse of notation we drop the tilde. In this limit the full S-matrix and consequently all 
S-matrices entering the Bethe-Yang equations become of difference form as appropriate for 
a relativistic theory. 

At the level of the simplified TBA equations this relativistic limit is implemented 
in a very simple fashion. In the interpolating theory have the following three types of 
convolutions 

/oo 
dtf(u,t)s(t-v), (5.1) 
-oo 

dtf(u,t)s(t-v), (5.2) 

-u b 

/-U b POO 
dtf(u,t)s(t-v)+ dtf(u,t)s(t-v). (5.3) 
-oo J u b 
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where 

s{u) = i^bf ' (5 - 4) 

and 

Ub = — log = — arcsmn g sin — . (5.5) 

irg 1 — t, ng \ kJ 

Rescaling the rapidities and taking the infinite coupling limit the points ± g go to positive 
and negative infinity respectively, and consequently the limit can be summarized as 

f*s->f*s, (5.6) 
f*s^f*s, (5.7) 
fis^O, (5.8) 

where on the right hand side we of course have the properly rescaled kernel 

s <»> = 4dhr < 5 - 9 > 

As a quick example of this simple consideration, let us consider the mirror energy satisfying 
the following identity 

£ M (K + 1)^ = Sq,J*s . (5.10) 
This implies that in the relativistic limit we should have 

£i-£ 2 *s = 0, (5.11) 

where in this limit 

lim 5 ^)=2cosh^^S. (5.12) 
9-s.oo \9 J 2k sin r 

Since Eq is a meromorphic function on the whole n-plane in this limit, the identity above 
nicely follows from 

£t+£i=£2- (5.13) 

In this limit the discontinuities of the Y-functions following from the TBA equations 
have disappeared completely. The dressing phase in particular only enters in defining 
the asymptotics of the Y-functions but does not enter in the TBA equations explicitly 
anymore. This is because before taking the relativistic limit the dressing phase enters 
the simplified TBA equations with a * convolution, as we have carefully proven for the 
deformed dressing phase in appendix D, which goes to zero in the relativistic limit. The 
only 'coupling constant' left in the game is the level k, and as in the interpolating theory its 
effect comes in through the boundaries on the Y-system itself rather than the discontinuity 
relations the Y-system is supplemented with. It should not be surprising that the kernels 
in the simplified TBA equations do not depend explicitly on the coupling constant; this 
is a common feature of relativistic models. The information on the value of the coupling 
constant comes in through the kernels of the canonical TBA equations which do depend 
on it, which in turn define the asymptotics of the Y-functions. 
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There is one special case however, which is Y + . Before taking the limit, Y + can be 
obtained as the analytic continuation of Y— through its cut along the real line. Now in 
the relativistic limit this cut disappears and Y + is a completely independent function that 
truly does not have a Y-system. It can be expressed in terms of Y_ and Yq as 

Y + {u) = Y_{u)e l < 1+Y Q> K Qy {u) , (5.14) 
where in the relativistic limit 

Id sink ^(u-iQ) cosh ^(u + iQ) 

K Qv( u ) = ir-.-r lp g c - v, *i ■ n \ — vT^w vtY • ( 5 - 15 

zm au sinh j^{u + iQ) cosh (u — &<yj 

In general the relation between Y + and Y_ will have state dependent driving terms because 
of the convolution involving Yq. Applying s~ l to the equation for Y+ does not appear to 
give more insight. 

Apart from this subtlety, the Y-system becomes a finite set of algebraic relations 
between meromorphic functions on the u-plane. Of course, also in this limit the Y-system 
depends on the state and twist in the manner described above. 

Transfer matrix. To obtain the tranfer matrix in the relativistic limit we rescale the 
rapidities u — >■ u/g, — Because of this, the factors Q 11 and Q m are basically unchanged 

K A 

Qi^ J] sinh ^(n - + Bi) = Q%. (5.16) 

i=l 

Similarly, let us introduce the following functions for the K l roots 

K 1 K l 
Q l B = Y[sinh^-(u- Ui + Bi), Qb = JJcosh^-(u -u { + Bi). (5.17) 

i=i i=i 

Then, we find that, after discarding the overall factor n^=i \j f^' ^ e transfer 

matrix (3.22) reduces to 

. \ " itt(a-2m) ^«-2m+1^2m-l Wt a \^ ia(a-2m) ^«-2m-l ^a-2m+l ^-a 

m=0 ^a-2m m= i ^ a+ i^ a _i ^a-2m 

(5.18) 



e ia(o-2m-l) ^tt-2m-l^a-2m-l 



m=0 



o n o m o n o m 

^a-2m+l , -0 ^a-2m-2 

•*-a— 2m ^a— 2m— 1 -^a- 2m— 2 -^a— 2m— 1 



Analogously, we derive the following expression for (3.30) in the relativistic limit 



e 



s 2i/3mQU glU Ql s i P{ 2m-l)Qll gill qIU 

2^ ~QUI qui e ~q~1 Z-> QUI QUI qui h 

m= "^s— 2m— 1 -^s— 2m+l -^s— 1 m= i -^s— 2m— 1 -^-s— 2m+l 

(5.19) 



Ql Q 1 1 ^ e 2 ^ m Q n 9 QP O™ . Q 1 . ^ e ^(2m+l)gH _ Q] u 



m=1 ^s-2m-l^-2m+l^s+l «-s+l m=0 ^s-2m-l^ 



s— 2m+l 
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For completeness, let us also give the special functions t used in the asymptotic solution 
Q} ,Q l , Q 11 



, e ia(a+2) Q-a-lQ-a+l Qq 



(5.20) 



*(a+l) Q-a-lQ-a-1 



o u o m o n o in 



Qi-,0 1 ,Q m ,Q n i , 



a-*-— a— 1 a— 2 a— 1 

q\ Q ii Q m m 



+ 



o 1 o 1 o 11 o 111 o 111 o 111 o 1 o 



(5.21) 



in 

s+l 



y s _iy_ s+ i wi_ s wi_ s _3Wi_ s _iWi s+ i 

Note that the transfer matrices are of purely trigonometric form. 



s+l 



6 Concluding remarks 

In this paper we discussed the effects of excited states on the quantum deformed AdSs x S 5 
mirror TBA equations. Doing so, we uncovered an interesting feature whereby the Y- 
system depends on (the excitation numbers of) the state under consideration, albeit in a 
mild way. This feature depends crucially on the root of unity deformation and the fact that 
we have a nested system, and to our knowledge has not been observed before. Similarly 
our TBA equations and Y-system depend explicitly on twisted boundary conditions, but 
this dependence is not unlike that in the XXZ spin chain, replacing twists by a magnetic 
field. 

In the asymptotic limit these interesting results can be verified through our asymptotic 
solution given in terms of the transfer matrix of the ^-deformed Hubbard chain and the q- 
deformed mirror bound state dressing phase, both of which we have explicitly constructed. 
This asymptotic solution also allows us to construct excited state TBA equations. It would 
be interesting to see to what extent the deformation can qualitatively affect the analytic 
properties of the Y-functions for an excited state in addition to the new effects which remain 
at the level of the Y-system as discussed above. As we saw in this paper for example, the 
ground state Y-functions of the twisted deformed model are nearly identical to those of the 
undeformed twisted model. 

Our equations reduce to a particularly simple form in the relativistic limit, where they 
are conjectured to describe the Pohlmeyer reduced superstring theory. As these equations 
still have a coupling constant dependence in the form of k, this limit is nontrivial and would 
be interesting to study in its own right. At this point however, it would be instructive to 
first understand the relation between the ^-deformed S-matrix and the (perturbative) S- 
matrix of the Pohlmeyer reduced superstring in more detail. Interestingly, there might be 
various subtleties in the relation between these S-matrices to which the resulting Bethe 
and subsequent TBA equations would be (rather) insensitive. 

From the point of view of finite size integrability it would also be interesting to consider 
the case where the deformation parameter q is taken to be real, where the theory shows 
some interesting features already mentioned in our previous paper. From the point of view 
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of condensed matter physics, both types of deformations are interesting deformations of 
the Hubbard and related models. We can address the thermodynamics of these models in 
the spirit of e.g. [71] with minor modifications of our equations. 
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A Representations and S-matrix 

The asymptotic solution of the Y-system will be formulated in terms of transfer matrices 
which partly take value in bound state representations. In this section we will set some 
notation and discuss these representations together with the corresponding S-matrix. 

A.l Bound state representations 

The quantum deformed S-matrix respects the quantum affine centrally extended psu g (2|2) 
algebra. This algebra is generated by four sets of Chevalley generators Ei, Fi, Ki(= q Hi ) 
and two sets of central elements Uk, Vu (k = 2, 4) with Uk being responsible for the braiding 
of the coproduct. The generators 1 and 3 form two quantum deformed su(2) algebras. 

The symmetric matrix DA and the normalization matrix D associated to the Cartan 
matrix A are 



DA 



( 2 -1 -l\ 
-10 10 
1-21 
\-l 10/ 



£> = diag(l, -1,-1,-1). 



(A.l) 



The algebra is then defined by the following non-trivial commutation relations, 



KiEj 



q DA *EjKi, 



KiFj = q~ DAij FjK{, 
{E 2 , F A } = -~g6T x (K 4 - U 2 U 4 1 K^ 1 ), {E 4 , F 2 } = ~ga +l {K 2 - UdJ^K^ 1 ) 



[E j ,F j } = D 



33' 



[E i ,F j } = 0, i^j,i + j^6. 



(A.2) 



These are supplemented by Serre relations (j = 1, 3, k = 2, 4) 

[Ej, [Ej,E k }} - (q - 2 + q-^EjEkEj = 0, [E X ,E 3 ] = E 2 E 2 = E 4 E 4 = {E 2 , E 4 } = 0, 

[Fj, [Fj,F k ]} -(q-2 + q-^FjFkFj = 0, [F 1} F 3 ] = F 2 F 2 = F 4 F 4 = {F 2 , F 4 } = 0. 

(A.3) 
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The central elements are related to the quartic Serre relations (for k = 2, 4), 

{[Et, E k \, [E 3 , E k }} - (q - 2 + q- 1 )E k E 1 E 3 E k = g k a k (l - V k 2 U 2 k ), 

{[F lf F k ], [F 3 , F k }} - (q - 2 + q- 1 )F k F 1 F 3 F k = g k a k \V k 2 - U k 2 ). (A.4) 

In total, this algebra has three central charges, 

d = K^lKs, 

C 2 = {[E2,Ei], [E 2 ,E 3 ]} -(q-2 + q~ 1 )E 2 E 1 E 3 E 2 , 

C 3 = {[F 2 ,F 1 ],[F 2 ,F 3 )}-(q-2 + q- 1 )F 2 F 1 F 3 F2. (A.5) 



Finally, the central elements V k are defined by the relation K l K k = V k . 

Hopf algebra. The coproduct of the group- like elements X 6 {l,Kj,U k ,V k } (j = 
1,2,3,4 and k = 2,4) is defined in the usual way, A(X) = X (g) X. The remaining 
Chevalley generators have a coproduct that is braided by the central elements U k 

A{Ej) =Ej®l + K^U^U^ 1 ' 4 8) Ej , A(F j ) = Fj ® Kj + U~ 5j ' 2 U^ 4 (8) . (A.6) 

Representation. Following [64] will use the g-oscillator representation to describe the 
bound state representations (short symmetric represenatations) . The bound-state repre- 
sentation is defined on vectors 

\m,n,k,l) = (4r(alr(a\) k (4) l \0), (A.7) 

where the indices 1, 2 denote bosonic and 3, 4 - fermionic oscillators. The total number 
of excitations k + I + m + n = Mis the bound-state number and the dimension of the 
representation is dim = 4M. This representation constrains the central elements via U := 
U 2 = U^ 1 and V := V 2 = V^ 1 and describes an excitation with quasi-momentum p related 
to the deformation parameter as U 2 = e ip . 

The explicit action of the triples corresponding to the bosonic and fermionic sl q (2) in 
this representation are given by 

Hi\m, n, k, I) = (I — k)\m, n, k, I), H 3 \m, n, k, I) = (n — m)\m, n, k, I), 

Ei\m, n, k, I) = [k] q \m, n, k — 1, 1 + 1), E 3 \m, n, k,l) = \m + 1, n — 1, k, I), 
Fi\m, n, k, I) = [l] g \m, n, k + 1, 1 — 1), F 3 \m,n,k,l) = \m—l,n+l,k,l). (A. 8) 

The supercharges act on basis states as 

H 2 \m,n,k,l) = — /c l|m,n,A;,Z}, 

E 2 \m, n, k,l) = a (— l) m [/] g \m, n + 1, k, I — 1) + b \m — 1, n, k + 1, 1), 
F 2 \m,n,k,l) = c [k] q \m + 1, n, k - 1, 1) + d (-l) m \m,n- l,k,l + l). (A.9) 

Here [n] q = (q n — q~ n ) /{q — q~ l ) denotes the g-number and C is the g-factor of the central 
element V = q c and represents the energy of the state. 
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The explicit parametrization of the representation labels in terms of the conventional 
x^ variables is 



g g ax 

-7, b 



[M\ q " \j[M] q l x- ' 

M M _ _ 

_g 7 igq 2 d= I g gg^Vx^-x (aw) 

[M] q aV g{x+ + i) ' y [M] q i 91 £x+ + 1 ' 1 ' J 



The central elements in this parametrization read as 

g M x- + $ q x~ ix+ + 1 ' g A/ Car - + 1 9 x" x+ + £' 



while the shortening condition becomes 

^(^ + ? + f + i)-«"(*- + F + < + ?)' (A ' 12) 

here £ = —ig(q — q' 1 ) and g 2 = g 2 / (1 — <? 2 ((7 — f?" 1 ) 2 ). 

The action of the affine charges H4, E4, F4 is defined exactly the same as the regular 
supercharges, but is subject to the following substitutions C — > —C and (a, b, c, d) — > 
(a, b, c, d). The affine labels a, b, c, d are acquired from (A. 10) by replacing 

V-^V' 1 , x ± ^^r, 7^^?, a-^aa 2 . (A.13) 
x ± x + 

Let us introduce the multiplicative spectral (evaluation) parameter u and a multiplicative 
parameter u of the algebra 



1 - U 2 V 2 

v 2 ~u 2 



girit 



ek. (A.14) 



To keep the the expressions in the remainder of the paper manageable, we introduce the 
following short-hand notation for the basis vectors of the representation 

\e k ) = \0,0,M-k,k), |e 3 , fc ) = |0, 1, M - k - 1, k), 

\e Ajk ) = \l,0,M-k-l,k), |e 34 ,fc> = \1,1,M -k-l,k- 1) (A.15) 

and we suppress the total bound state number M. 

A. 2 S-matrix 

We consider the bound state S-matrix which is an intertwining matrix of the tensor space 
spanned by the vectors 

\mi,ni,ki,h) ® \1ri2, n 2 , &2, h). (A. 16) 

Here < mi, ni, rri2, ri2 < 1 and ki,li,k2,h > denote the numbers of fermionic and 
bosonic excitations respectively with the bound state number Mj being the total number 
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of excitations M, = rrii + + ki + l{. The S-matrix is required to be invariant under the 
coproducts of the affine algebra, 

SA(J) = A op (J)§, for any JeQ. (A.17) 

The invariance under the bosonic symmetries AH\ and A.H3 requires the total number of 
fermions and the total number of fermions of one type 16 

Nf = mi + m 2 + m + n 2 + 2l\ + 2/ 2 , 

JV /3 =mi+m 2 + HI 2 . (A.18) 

to be conserved. This conservation divides the space (A. 16) into invariant subspaces de- 
noted I, II, III in [64]. The explicit expressions for the coefficients are given in [64]. In this 
paper we will restrict to the case where M 2 = 1 and we have for completeness included the 
S-matrix elements that are needed in our derivations in Appendix B.l 



B Algebraic Bethe Ansatz 

In this appendix we give various details concerning the algebraic Bethe Ansatz preformed 
in Section 3. 

B.l Elements of the S-matrix 

In order to derive the transfer matrix we need the following elements of the S-matrix 
where the first leg is an M-particle bound state and the second leg is in a fundamental 
representation. The S-matrix that scatters subspace I is 



5 Note that a bosonic excitation is interpreted as a combined excitation of two fermions of different type. 
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For subspace II we have 



Finally 



-m;l 



X 



m,0;l 
m:3 



'm;3 



'm;6 



X 



m,0;3 
m;6 



5! 



m,0;l 
m;6 



'm;l 



m,0;6 
m:3 



^m;l — 



^m,0;l 



m 2 x + 



x+ q M ~ 2m m 



<?u 2 



^2^2 X x — X 2 

[M - m] q 7i x\ 



gu2 



^m,0;2 



^m,0;2 
^m;2 



y mA 



OymflA 



^m,0;4 
^mA 



72 »i 



[M] g 7i £+ - x 2 
m x7 — Ixi 

q 2 -r : 



X i 37 



1 ^2 



q m a 




\ x l x l )\ x 2 x 2 )\ x l x 2 Ml 



7i7 2 7[M] 




1-9 



[m] q [M 



M-m 



m 



[m] q (xf - x 1 )(x 1 x^x+ - x 2 ) 

[M]q (xf - X 2 )(x^xfx 2 - Xi) 
„M 



\q ^2 ^2 



1 37-^ ^ (*^2 ^ 

[M],(x+-x 2 -)(1-x+x+Kt| 



x7x 



_(x+ -x 2 )(l 



M+l X x — X 2 1 



X 



,+ 



x 2 1 




^M-m [ m ]<? x l )( x l x j x 2 X 2 

[M] q (xf - x 2 )(x ] ;x|x 2 f - xj~ 



x 1 x 2 



X -i 37 



37 1 37 n 



2 i x x x 2 



l 7! [M — m] g X 2 — X 

72 ^[M], x+ - x 2 1 - x^x-j 
a (x^ — x 2 )(x^ — x^)(x 



7pWTg7i72 (x+-x 2 )(l-x+x+) 




72 



37-i 37r 



^/\M] q li (x+ - x 2 )(1 - x+x+ 



3M 

g 2 



m M?7i72 ^ 



\/[M]q a x\ - x 2 \ - x\x 




x 1 x 2 



(B.2) 
(B.3) 

(B.4) 

(B.5) 
(B.6) 

(B.7) 

(Bi 



(B.9) 
(B.10) 
(B.ll) 
(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 

(B.17) 
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B.2 Commutation relations 

Here we list the needed commutation relations. Greek indices run over the fermionic indices 
3,4. The diagonal elements Tj depend on the spectral parameter z and the representation 
is described by the parameters without indices. The creation operator depends on the 
spectral parameter A and we further denote x + (A) = y and v = qX 



«(A) 
«(A) 



cc + z£SM(A) 



pP,mrpl 



aV — X 

! 2 : 

y — ar 



ay — x 
y — x i 



ay — x 

'■ 2 j 

y — ar 



ay — x 

n : 

y — 2P 




r c [(A)r + ... 




u — V. 




q u 



ri(A)CC + T3 3 4 4 ; 



34,mx 



+ 



(B.18) 
(B.19) 
(B.20) 




u — v 



'u — g'"v 



„7<5 /_M-2m-l 



u,A)Ti(A)T^, 



(B.21) 



where the dots stand for additional terms that do not affect the eigenvalue but should 
vanish upon solving the auxiliary Bethe equations. The auxiliary S-matrix r is given by 
(3.20). 

B.3 Auxiliary Bethe Ansatz 

Consider a general £)t(l) 2 invariant S-matrix 

fn(x,y) \ 

r 2 (x,y) r 5 (x,y) 
r G (x,y) r 3 (x,y) 
V n(x,y)J 



r(x,y) 



(B.22) 



satisfying the Yang-Baxter equation and unitarity (i.e. ri2(^i, ^2)'"2l(^2) ^l) = !)• Fol- 
lowing the algebraic Bethe ansatz approach we introduce the monodromy (M), transfer 
matrix (T) and vacuum state (|0)) on a chain with L sites 



M 



A B 
C D 



T = ttM = A + D, 



|0> 



8=1 



o, 



(B.23) 



The action of of the different elements of the monodromy matrix on the vacuum are readily 
computed 



A(q)\0)=l[r 1 (q,p l )\0), D(q)\0) = f[ r 3 (q, Pl )\0), C|0) = 0. 



(B.24) 
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And the operators from the monodromy matrix are found to satisfy the following commu- 
tation relations 

A(x,m = W) B{v)A(x) - ^) B{x)A{y) ' (R25) 

D(X)B{V) = r ^i D(y)B[x) - r ^M, D(x>B( « l (R26) 

B(x)B(y) = r -lpHl B (y)B(x). (B.27) 
r 4\x, y) 

For a state with K excitations we make the ansatz 

fl(Ai)...B(A*)|0). (B.28) 

Then putting all the ingredients together we find the following eigenvalue for the transfer 
matrix 

a(«x) - n«n|f+M^ ^ 

supplemented with the Bethe equations 

L K 

1= r 5 (\ k ,q)r 3 (q, A fc ) _ -r-r ri(\ k ,Pi) rr n(Xi, A fc )r 3 (A fc , Aj) (B 30) 

^3(Afc,g)r 6 (g, A fc ) ^r3(A & ,p i )^r 3 (A i ,A fe )r4(A fc ,A i )' 

where we point out that the LHS is actually constant due to unitarity. 
C Dualization 

There is a way to relate the transfer matrices and Bethe equations in different gradings 
that goes under the name of dualization. They are related to each other by a change in 
variables and quantum numbers. Here we will apply a duality transformation to the y 
roots, giving transfer matrix in the su(2) and sl(2) gradings. 

The transfer matrix in the sl(2) is given by (3.22). Let us introduce the polynomial of 
degree K l + 2K 111 in terms of the multiplicative evaluation parameter 

K 1 _ + K™ l K 1 _ _ K UI 

p{ y ) = ns=^ iiffa - ^ - ^ ns=^ n^w* - ^ = pm - p^ (d) 



i=l i/I- i=l i=l i/Xj i=l 

where v = y+y ^_f^^ — . By construction the polynomial P has K u roots corresponding 
to the roots of the auxiliary Bethe equations (3.23). Factorizing P then yields 

P(y) = cl\(y-y l )Yl(y-y i ), (C.2) 

i=l i=l 
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where are the dual roots and there are A 11 of them. As a consequence we find the 
following identity 

K 11 K u 

n yi-a _ P{a) jj jjj-b 

Dualizing means to write the transfer matrix and Bethe eqautions in terms of y and A . 
For conciseness we will present he deriavtion in terms of the multiplicative parameters. 

Auxiliary Bethe equations. Let us now dualize the Bethe equations. Since y is a root 
of P it satisfies the same auxiliary Bethe equations (3.23) 

ifi ~ + r^K™ 

1 = rrJ_ y*-*i An ^ , (c.4) 

f = \ Vvm-Xi V X T 7=1 v Wi ~ Yk 

Next, we continue by looking at the dualization of (3.24). Introducing two auxiliary pa- 
rameters LU^ 

i ^+ + ^ + g + f + ^ + e + | 

w* = ~ = ~ q ^1 , (C5) 

and using (C.3) then allows us to write 

A ' n RU 1 (Vi ~ vtKVi ~ \) 1 P(ut)P(^) Rn 1 

l\ - gv< l\ q ( W - a£ )( W - ^" ^K")^(^) ? - V 1 ' ' 

From the definition of P it is straightforward to prove 

resulting in the following dualized version of (3.24) 

gw fc - Vj w fc - g^Wj 

where we used that A 11 = A" 1 - A 11 + 2 A 111 . From this we see that the dualized auxiliary 
equations can be simply obtained from (3.23) and (3.24) by replacing y -H- y and A 11 A 11 
respectively. 

Transfer matrix. After discarding the overall scalar factor and rearranging some terms 
in the transfer matrix, we find 

^ - - p^) + i^r y ^ + ( } 



i=l * - J ^ x- 

K 1 i 1 K 11 K™ 



p _ a-1 p/ T [a-2m]\pC 1 \ K l a _ m „m+l„ # K „„a-2m+l 

-K a ^ ^ J- r l 2 .[a-2m] J T-T(? U~ (? Uj-pT (/ U — Vj TT 11 '? 

^ y T3, ( ~.\a—2m\\T>. ( 1 '\ 11 rf,, _ ml . 11 n a-k-,, _ „m„. 11 



w,: 
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We can write this in a more convenient form in terms of the dual quantum numbers by 
using the following properties of P 



K 1 a -2m.. . D I [a- 2m] \p ( 1 \ X m , a+1 x2 

rr g u — gu t _ -rnx /^U^a-H ^ tt / q u — w t \ 

11 ^U-^Ui p i ( x +)p i (0 F ) 11 VgO-2m+l u _ w J ' V • J 

S g Q u-v, _ 1 ^ + )^) S g a - 2m u-v, 

11 q a-m u _ q m v . m (K l +2K lll -K 11 ) p( x [a-2m])p( 1 ) 11 g a u _ y . ' I " J 



leading finally to 

P(x+) P 2 (x-) 



t a 



P 2 (x+) P(x-) 



P(x-) P(^) yr u - q a Vi yr (g a - 1 u-w t )(g a+1 u-w t ) 
P2(x~) P!(Jp) 11 g a u - v< 1 = 1 (u - g^w^u - q a+l Wi) 

(C.12) 



P(^) P(x-) ^4 ^ g"~ m u - g m v, yj- 1 (g Q+1 u - wf)(g^u - w 



W; 



Pi(^+) P2(x~) ^ ^ g a u - v, 11 (g«~™u - g-+ 1 w i )(g a - m u - g' 
This is the dualized transfer matrix. 

Relating t aj \ to t\ )S . Let us introduce the polynomial P, which is the complex conjugate 
of P 

K 1 _ A 111 A" 1 + K UI 

p(y) = ii* (*< - x ) - 4^ n*- - v ) = ^ - ( c - 13 ) 

i=l Jx i i=i q q 2 i=i Wx+ i=i 9 

By similar arguments it can be shown that T\ s can be written as 

P(x+) P(^) ^r g^u-Vj yr u-q s+1 Wj gu - q s Wj 
1,3 P 2 {x+) P^) l = lu-g s Vi 11 q s u-qwi g»+lu - w< 1 ' J 

P(^)P(X+) ^^ g^u-g^ ^j- 1 ( U -^+l Wi )( U -g^l Wt ) 

Pl(^)P2(ic + ) M u-g s Vj 11 (g s - m u-g m+1 w i )(g s - m u-g m - 1 w i ) 

771 — 1 I — _1_ I — 1 

It is now trivial to see that t\ tS maps into £ a> i (after changing k o s — k) under complex 
conjugation (i.e. map (a; , g) — > (x^, 1/g)) up to an overall normalization. 

D The dressing phase 

D.l The dressing phase for fundamental particles of the 'string' theory 

The crossing equation that follows from the g-deformed i?-matrix has a solution which is a 
natural deformation of the AdSs x S 5 dressing phase [65-67]. This solution, which we will 
denote a, was found in [44]. The dressing phase in our conventions, a, is related to a as 

2/ \ -2/ \ x f x 2 x i £ x t ~^ £ _ ~2( \ P i x l) /t-> i \ 

cr (xi,s 2 ) = o- (^1,^2 — —n , — _ — = cr ^1^2)^7 — r • D.l 

xj x^ + £ x 2 + £ P(x 2 ) 
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Both a and a solve the AdSs x S 5 crossing equation [68] in the limit q — > 1. The dressing 
phase a is conventionally written in the form 

a(zi,zi) = e^ (2l ' 2l) =expi( X (xt,x+)-x(xi,xt)-x(xf,X2)+x(xi,X2)) > ( D - 2 ) 
where when both particles are in the string region, the x functions are given by 

dz I J 1 duu , r ? 2(l + y(u(z) — u(w))) 



y(x 1 ,x 2 ) = i <p <p log^-^ 2 (D.3) 

7\z\=\ 2m z - x x J H=1 w - x 2 2m & T q2 (l - f(u(z) - u(w))) 

Under analytic continuation to other regions of the rapidity torus the expression for the x 
function changes. The above double integral is commonly denoted 

dz I £ 1 dw , r g 2(l + ^(u(z) — u(w))) 



J\ z]=x 2mz-x l J ]M=1 w-x 2 2m B r, 2 (l-f(ti(z)-«H))' 

and is equal to the x function in the string region. The function has a discontinuity on 
the edge of the string region (\x \ = 1), and to properly define the x function beyond it 
we will need further terms, which leads us to introduce the function 

J\ z \=i 2m z - x 2 r ?2 (l - f(«i - u{z))) 

With these definitions we are ready to give the expression for the dressing phase in the 
currently relevant regions of the torus, defined as 

K : |x ± | > 1 , Tl x : < 1, > 1, K 2 : \x + \ < 1, |aT| < 1 . (D.6) 

7l a ,b denotes regions on the product of two rapidity tori in the obvious fashion. We will be 
most interested in the (bound state) dressing phase on the real line of the mirror theory, 
which lies in region Tin. For completeness let us first briefly repeat the explicit proof that 
a satisfies the crossing equation, given in [44]. In order to do so, we need to analytically 
continue the dressing phase to region 7^2,0- 

D.2 Proof of crossing 

The dressing phase in Tlip 

Continuing <I> through \x^\ = 1 to \x^\ < 1 gives 

Ki, : =$(x+,xf)-y(x+,xf), (D.7) 

x(xT,x 2 ) =$(xi ,xf). (D.8) 

The dressing phase in 7£ 2 ,o 
With \xf\ < 1 we have 



Tl 2fl : X (xf,xt)=Hx+,xf)-^(x+,xt) + ^log^ ±, (D.9) 

X(xr. x t) = HxT,xf) - , xf) . (D.10) 



1 






x 2 






x- 


— x 2 
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Identities I 

In order to prove the crossing relation between 7£ 2j o an d 7£o,o we w ih need some identities. 
First of all we have 1 

*(X1,X2) + *(1/X1,X2) =*(0,x 2 ), (D.ll) 
x 2 ) + *(a:i, 1/^2) = *(a;i, 0) . (D.12) 

Secondly we need 

^(l/x^,x^)-^(l/xf,x^) + ty(l/xf,X2) - $>(l/x±,x 2 ) = 
1 1 1 L_ 1 + J_ 

J 1 x ~ x + x + x + x + 



- log - — ^ — + f lo s — r • 

x 1 x 2 x 1 x 2 x 2 

This identity follows from the identity 



x 2 - X2 - -4- 

*(l/xi,x 2 )-*(l/xJ,x 2 ) = -ilog ^ (D.13) 

£ 2 ^2 + ? 

whose derivation plays an important role in the fusion of the mirror dressing phase below 
as well, so let us discuss it in some detail. To prove this identity, we begin by combining 
the two integrals 

^(l/x^,x 2 )-^(l/x+,x 2 ) = (D.14) 

= il ^->—\ ryq + |K - "(*))) r>(i - ¥« - "(*))) 

1 f [g[=1 2ni z-x 2 ° g r, 2 (l - f («i - «(*))) I>(1 + f («+ - u(z))) 

Next we use the defining property of the r g function to find 

f dz 1 1 — a^( u r-«0)) 1 — ( f90( z )-'«i~) 

,x 2 )-^(l/x+ x 2 ) = i <b — log 

J\ z \ =l 2mz-x 2 1 - q 2 1 - q 2 

where we note that the log(l — g 2 ) terms do not contribute to the integral as the integration 
contour can be shrunk to nothing. After integrating by parts 18 we find 

¥(l/a?i ,2:2) - ^(l/xf,x 2 ) = 

\og(z - x 2 ) (-1 + + + + _ 1 _ 1 ) 

1*1=1 27T« &V ^ * *-*f Z ~~^x 1 J 

X 2 - 4fr X 2 - 4r 

= -ilog ^ (D.16) 

£ 2 £ 2 + £ 

17 These identities follow by changing variables from z to z _1 in the second integral on the left hand side. 
18 Note again that |a?3 1 > 1. 
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by summing the poles within the unit circle. Here we used the fact that \xf\ > 1. Had it 
been different, the respective x would have to be replaced by l/x in the above expression. 
Note that ^ is invariant under inversion of its first argument since u(x) = u(l/x), so that 
the above also directly applies to 9{xJ,X2) — ty(xf,X2). 

In the construction of the mirror bound state dressing phase we will need a further 
identity 19 , namely 

*(l/xr,0) -fr(l/z+0) =tlog 2 +ilogxrx+ _ Xl -— l — (D.17) 

(i - <rr x 1 +iix x +1 

= i log J + % log x^xf - i log q-^-i/a) (1 - e 2 ) . 

which follows by carefully integrating (D.15), noting that now the constant log terms in 
the integral cannot be dropped. In the limit q — > 1, get the undeformed result 

a 2 

V(l/xi,0) -V(l/xt,0) =i\og~+i\ogxJxt (D.18) 

— f 2 a 2 

since for q ->• 1, £ ->■ while -)■ ^. 

Crossing 

Putting the above together we can directly prove crossing between 7£o,o an d 7^2,0 ■ This 
firstly makes use of the identity 

*(xi,x 2 ) + ^(l/xi.xa) = $(0, x 2 ) , (D.19) 

valid for \x\\ 7^ 1. Using this we find that 

A§ = 0(z 1 ,z 2 )+e(z 1 + uj2,Z2), (D.20) 

is given by 

AO =*(l/xj",a; + ) - *(l/xf + , x 2 ) - #(l/a;i 

+-iog x 1 r ~ x i^" X2 „ . (D.21) 

* x 2 x l x 2 

Then we use the identity 

y{l/Xi,X+)-V(l/x+,X+) + X2 ) - = 

1 l-^rl-^r 1 + 4 

-rlog - Xi ; 2 - ^ +ibg — x f. 

I 1 =^ 1 - -qr— =■ 1 + 4r 

19 In these derivations we will not be very careful about factors of iiv - the typical ambiguity in defining 
\ log(— l) 2 - as they do not affect logo -2 . 
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to finally find the crossing equation 



1 _ J i L_ _ + J__ X - i + _L 

1 . x, xi x t x t x l — X 2 X l ^ -1 X 2~ 



Afl = T log- ^ + T i og ^ ^ 3 + ^log 

1 _ _ 1 _|__ ^ _ Xo X~\ Xry 

1 



1 + 



log 



%±^h(x 1 ,x 2 ) . (D.22) 



In other words 



o-(xi,X2)a-(l/xi,x 2 ) = + ; ^ h(xi,x 2 ) ■ (D.23) 
If we now rewrite this in terms of <r we get 20 



^2 + 
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<r(ai,a;2)cr(l/ai,a;2) = --rrM^i) x z) > (D.24) 

which is of course precisely the crossing equation a is supposed to solve. 

D.3 The improved bound state dressing phase of the mirror theory 

We would like to construct the dressing phase for bound states of the mirror theory. In 
order to do so we take the approach taken in [69] in the undeformed case, and take the 
constituents of the mirror bound state such that only the first particle lies in region 1Z±, 
i.e. \xi\ > 1 but \xf\ < 1, while all other particles lie in region K with |x ± | > 1 [27]. In 
order to sum up all contributions to the bound state dressing phase of the mirror theory, 
in addition to the above we need the dressing phase in TZ\^. 

The dressing phase in T^-i i 

Continuing the x functions to the region with \x±\ < 1, \x±\ > 1, \x% \ < 1 \ x %\ > 1 
means we have to add additional terms to the expressions in region T^o, corresponding to 
crossing = 1. This means the two <E» functions give extra ^ contributions, but also the 
already present \& function gives an additional contribution. The full x functions are then 
given by 

72-1,1 : x( x t, x t) = x t) + ^i x 2' x l) - ^( x f, x t) 

i>(i + !(*+ + ^-4-^)) 

+ i log ; 



i>(i -f(*+ + ^-4-^))' 

X( x t, x 2~) = $( x l> x 2) - ®( X 1 ) X 2 ) > 
X( x l, x t) = $( x l, x t) + %( x ti X l) ' 

X( x i, x 2 ) = $( x i, x 2 ) > (D.25) 



20 Note that the dispersion implies x + + — - = \ 



_ 4i _ 



The bound state dressing phase I 

As already indicated, the constituents of a mirror bound state can be taken to lie in 
the regions 1Zi,i, 7^1,0, 7^o,i ; and TZo,o- By antisymmetry of the dressing factor the above 
sections provide the appropriate dressing factors for the constituents, which can then be 
fused. At this point we would like to introduce the improved dressing phase £ 

X(x 1 ,x 2 )=a(x 1 ,x 2 )- ^ (D.26) 

1 _ 4- 



1 i_ 

1 — = + 



where P is defined in (D.l). The chosen representative constituents for a Q-particle bound 
are parametrized as 

xj(u) =x s (u + (Q-2j)-), j = l,...,Q, (D.27) 
J 9 

x+(u) = - 1 , xf(u) = x s (u + (Q - 2j + 2)-) . (D.28) 

x s (u + iQ/g) J g 

We then denote the final mirror theory bound state parameters as 

vt = x i/q( u ) > yf = x yQ>( v ) ■ ( D - 29 ) 

Summing up all contributions for the Q-particle bound state - Q'-particle bound state 
improved mirror dressing phase using the appropriate form of the \ functions in each 
region, we directly find 

-ilogS QQ '(yi,y 2 ) =$(yt,y£) -$(yt,V2 )-*(vT> V2) + Hv!>V2) 

- *(vt, vZ) + *(yf,y 2 ) + ^(yt,yt) - ®(yZ, vx) 

1 - Q-l 1 - -J-r Q'-l 1 

'^T^=t^i^t^T^t- lD ' :,( " 

3/, vt J=l x, y+ k=l y z. 



Vl y 2 J =i x j Vi ft=i Vi z k 

where we emphasize that 

\yf\<l, \yi\>l, |y 2 + |<l, |i/J|>1. (D.31) 

This result has an apparent dependence on the bound state constituents. To manifestly 
remove this dependence we need a few more identities. 
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Identities II 

Firstly we have 

„- _ „+ V2 ~ IF ^ - ^= V2 ~ ■£= 

*(^ )2 / 2 -)-*( y r,%) -iog^— \ n - J ' (d.32) 



2 



y 2 +U V y 2 / V V1V2J j = i V x i y i J 

which follows by applying the defining property of T q Q times in the derivation (D.13,D.15,D.16) 
and using the fact that \ < 1 while 1 2/^2 1 > With both yi and y 2 referring to particles 
in region 1Z\ this also immediately implies 

\Vi +U \ Vx J \ 2/i 2/2 / J=i V / 
Then we also have the bound-state analogue of (D.17) for a particle in region 1Z\ 



^(yjf , 0) - 0) = zQ log ^ 2 + Hog TT (x- ) 2 - »Q log g 

(1 - g^ y+ AJL J 

By applying the basic property (D.12) we can reexpress terms of the type ^ , (y^,y 2 f ) an d 
apply the above type of identities to find 

y+) - *te,yt) = Vfa, li) " , ^) + *(Vi, 0) - *(l/r. 0) 

=^iogfe + + i)(i - 9 2 ) 2 (i - r a )ff"* (tt "^+ 

Q-l 



i log(yJ - 2/1 ) ( - ^ JJ {xj - yp 2 



^iogfe + + i)(i - g 2 ) 2 (i - r V* 1 * - ^ 



/i.>.M//. 2 -//r }(//.t -4) II fi-^r ) ( D - ;!;5 > 

f 1 / f=l V X j ^2 



1 



ilog JJ (r i _ £2 ( _ q-iaiv+Q'g) 

i=i 



Similarly we find 



21 



v(vt,yt) - Hv2,yt) =^'iog(^ 1 + + i)(i - g 2 ) 2 (i - r 2 )q iB(v ^+ 

y-i / . \ -2 



>iogw-v 2 -)(^-^)n (i-^j (D.34, 



1 Here the parameters x- of course refer to bound state number Q' and center v. 
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1 



Putting the above identities together with (D.32) we have 

Q-l 1 Q'-l 1 

-*(yt,yt) +*(yt,y2) + *(yt,yt)-*(yt,yl)-^n 1 II j "? fc 

= g i". vt) + % ) + *(va", yf) - *(yt: vl)) 
+ \ (-*(vl,yt) + *(fi~ 2/2" ) + *(i/J -f*) - *(% ,fD) 

Q x / J9(u+(Q _ 2i) i) _ -HKfO'J A /(1 _ 2) 

+ « log ) ; r-^ log ^ ^ 

Jl^-T 1 ( g _is( " +(Q ' _2i) s ) - j _ g2 ) 2 2/i 2/2 

+ ^-iogg 9 (C2/r + l ) — = ^-iogg 9 (a/2+1) — — 

Then we can use the defining property of the T q function to find 

■nO'-l ( -ig(v+(Q'-2j)-) -ig(u+Q-)\ o N / 

n?=i (g 3V ~g g3J j/(i-g 2 )r g2 (i-f(n(y + )-^ 2 + )) 

nji 1 ^-^«+«-^|) _ g-^^J)^ /(l _ g2) I> (l + I - u{yt)) 
nfc 1 g-^^ + W'- 2 ^^ I> (Q' - ^ („(y+) - u(y 2 +)) 



(D.36) 



n Q-l ^(^(Q-2,)-) r<?2 (q + iJL {u{y +) _ u( j,+)) 

I> (Q>- % (u(y+)-u(y+))) 
r ig(Q'-l)v q ig(Q-l)u q V* 2 V V^l ^ V«2 ^ 

r g2 (Q + f(n( yi + )-n(y + ))) ' 



The bound state dressing phase II 

By using the above identities, the improved dressing phase given in eq. (D.30) can be 
written as 

-ilogZW ( yi ,y 2 ) = <S>{y+,y+) - <S>(y+,y 2 ) - Qfo, y+) + y^ ) (D.37) 

- 2 (* (2/i + ' fa") + *(yr> y 2 + ) - *(fi", V2 ) - *(»r»% )) 

+ J 2/f) + , vt) - ^(yi, iff) - *(iff , iff )) 



(Q'-|( U (y+)- U ( y +)))l- 1 

i lor 



V2 fl % 



OT g2 (Q + I - u( % + ))) 1 " ^ V fl f2 



+ - i og gQ-Q'g-«KQ+Q'-2)(u-f) _ 

This version of the improved mirror bound state dressing phase manifests its proper fusion, 
as any apparent dependence on the bound state constituents has been removed. 
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To see that the improved dressing phase is unitary it is important to realize that 
the logT^ combinations entering in <3?, ^ and explicitly in the formula above, are not 
real. Rather their imaginary parts cancel precisely against the manifestly imaginary terms 
in the above prescription. One way to show this is to realize that the above expression is 
independent of the constituent particles, and so should agree with the improved bound state 
dressing phase obtained by fusing over a manifestly mirror theory conjugation symmetric 
bound state configuration. The latter expression can be easily proven to be unitary by 
using the crossing equation and the generalized unitarity of the dressing phase about the 
string line. We have also verified unitarity of the above expression explicitly numerically. 
The further relevant analytic properties of the dressing phase are covered in the section 
below. Note that this expression for the dressing phase manifestly reduces to the improved 
mirror dressing phase of the undeformed theory. 

D.4 Action of the discrete Laplace operator 

The mirror-mirror bound state dressing phase is a holomorphic function in the intersection 
of the region and the mirror region, which in particular contains the mirror line. This 
follows by the discontinuity relations for the functions in their first argument, which 
we discuss below. Holomophicity and fusion mean that the dressing phase is immediately 
annihilated by the discrete Laplace operator Ag/ p = Jg/ ps -1 — (<5g'+i,p + $q>-x,p) for 
P > 1, since we are far from any branch cut ambiguities of the arguments of the dressing 
phase. For P = 1 this is no longer the case; acting with the operator 

A Q , A (v) = 5 Q ,, 1 s- 1 (v)-5 Q ,, 2 , (D.38) 

only annihilates the mirror-mirror bound state dressing phase for v £ (— Uf,, Ub) 22 ■ This is 
the reason the Y-system for Y\ only exists on this interval of the real mirror line. Because we 
do not encounter any branch cuts of the arguments of the dressing phase for v G (— u^, it&), 
annihilation within this interval still follows immediately by holomorphicity. 

To show that the dressing phase is holomorphic in the intersection of and the 
mirror region, in addition to the formulae already derived we will need the discontinuity of 
the \l/ function ^f(xi, x%) in its first argument. The relevant discontinuity is at u\ = u+2i/g 
with u £ (—Ub, Ub). To find it we proceed as in the undeformed case [69] and integrate by 
parts to write the function as 

m( Xl ,x 2 ) = - 9 - I ^ : log( 2; - ;K2 )^^ g2 (l + |(« 1 - u (^)))+^ 2 (l-|( Ul -u(^)))) , 

1 

where ip q is the g-digamma function, the logarithmic derivative of the T q function. From 
the defining relation of the T q function it is easy to see that the q-digamma function 
still has simple poles at the negative integers with residue negative one. This makes the 
discontinuities immediately clear as two poles 23 hit the integration contour when u\ = 

22 The significance of the interval (—Ub,Ub) lies in the fact that x m (u) covers the unit circle precisely as 
u runs over this interval. 
23 Recall that u(z) = u(l/z). 
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u + 2i/g. The corresponding discontinuity is then simply 

lim *(e e xi, x 2 ) - y(e- e x 1 ,x 2 ) = -i log ^ff - * 2 , (D.39) 

following from the two pole contributions just inside and just outside the unit circle 24 
respectively. Analogous formulae apply for discontinuities at u\ = u + 2ni/g, exactly as in 
the undeformed case [69] . For example for the discontinuity through u\ = u — 2i/g we get 

lim^{e e x 1 ,x 2 ) - ^(e- e Xl ,x 2 ) = ilog Xm ^ ~ X2 , (D.40) 

x m (u) X2 

where the different sign arises because the poles cross the integration contour in the opposite 
direction. This shows that the cuts of combinations like ^(y^yf) + ^{y 2 ,yf) precisely 
cancel in the intersection of Ri t i and the mirror region, making the dressing phase a 
holomorphic function there as mentioned above. 

With this identity we can also directly compute the action of Aq/ } x(v) on TP® (u,v) 
almost identically to the undeformed case [70], and show explicitly that it indeed vanishes 
when v G (—u^^Ub)- Acting on the functions, noting that the \& function is precisely its 
discontinuity when its second argument crosses the unit circle (cf. eq. (D.4)), we get 

($(yt,vt) ~ HVi, V2 ) - HViiVt) + HVi,y 2 )) Aqm(w) = *(x m (v),yt)-V(x m (y),yf) . 

(DAI) 

For the first line of the ^ functions in (D.37) we get 

-\ [nvt'Vt) + *(yi,yt) - Hvt^) - Hvi^ )) Aqm(«) = (d.42) 

r2 fa -m u - v )) 



log q -ww-')(u-v) SQ + - log 



2 q(2 + f(u-, 

which simply follows from the integral representation for 'F of eq. (D.5) when its second 
argument crosses the unit circle. A effectively acts on the first argument of the second line 
of ^ functions, giving 

- (y{y+,yt) + y(y 2 ,y+) - *(yt,Vi) ~ .vD) &q>,i(v) = 

= *(xm(v),yt) ~ ifc) ~ i V i~ X i V)V l~^ ■ (D.43) 

2 Vl ~ H-Jy) Vl ~ X m(v) 

where we used the discontinuity of the function we just computed, given by eq. (D.39). 
This can be partly rewritten as 

* , yt-x m (v) Vi - ^jv) ., vt ~ ^t{v) \y{ i. -act \ (t^aa\ 

- ^og— — = i\og— —\ — p + ^logg w Sq(u,v), (D.44) 

2 w - ^jv) yi - Xm ^ yi - ^) V 2 



'Effectively giving x m and \jx m respectively. 
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which follows from the general identity 

\~^ X i~ X i = q - {Q+M) S Q - M , (D.45) 
1 +Sr xj - xj 

where particles i and j have bound state numbers Q and M respectively. Finally the last 
terms give 



(Q'-IK^)-^)))!-^ 



iQ'T q2 (Q + f (u(yt) - u(y+))) 1 " ^+ V Vi W 
^2 / Q _ »9 

^ 2 2 



-i log J- - i log _ x i v \ % (D.46) 



and 




L - logqQ-Q'q-^+Q'-^-^AQ,^) = - logq^q-^- 2 ^-^ (D.47) 
respectively. Adding everything up we get zero as promised. In summary, this means 

Z QQ 'A Q , A (u,v) = 0, for v E (—Ub, Ub) . (D.48) 

When v is outside this interval the result of applying A is nonzero. In the undeformed case 
the resulting kernel can be cast in a simple form [70], we will not pursue a simple version 
of the resulting deformed kernel here. 
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